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Dynamic nuclear polarization (DNP) is used in nuclear magnetic resonance (NMR) to transfer 
polarization from electron spins to nuclear spins. The resulting nuclear polarization enhancement 
can, in theory, be two or three orders of magnitude depending on the sample. In solid state systems, 
however, there are competing mechanisms of DNP, which, when occurring simultaneously, reduce 
the net polarization enhancement of the nuclear spin. We present a simple quantum description of 
DNP and apply optimal control theory (OCT) with an open quantum system framework to design 
pulses that select one DNP process and suppress the others. We demonstrate experimentally an 
order of magnitude improvement in the DNP enhancement using OCT pulses. 


I. INTRODUCTION AND MOTIVATION 

Dynamic nuclear polarization (DNP) has long been 
used in NMR to increase the signal to noise of an ex¬ 
periment by transferring polarization from electron spins 
to nuclear spin^. The thermal polarization of an NMR 
sample is proportional to the external magnetic field, Bq, 
the gyromagnetic ratio, 7, of the spins in the sample, 
and is inversely proportional to the temperature, T, in 
the high temperature limit, P « 'yBo/ksT. The electron 
gyromagnetic ratio, 7e, is two to three orders of magni¬ 
tude larger than 7„, the gyromagnetic ratio for a given 
nuclear spin. Thus electron spins are more highly polar¬ 
ized at thermal equilibrium than nuclear spins. DNP has 
become a valuable tool for NMR spectroscopy, medical 
imaging, magnetic sensors, and semiconductor spin stud¬ 
ies due to the increased sensitivity provided by nuclear 

polarization enhancement^2lHl_ 

While DNP has been studied extensiveljO^, a full 
quantum description has been missing from the litera¬ 
ture. In fact, DNP offers a illustrative example of open 
quantum system dynamics. Here we consider a two-spin 
electron-nuclear hyperfine coupled spin system, in which 
the polarization transfer can occur through several path¬ 
ways. Some of these DNP processes require decoherent 
mechanisms, so it is useful to consider the pathways of 
DNP to be various information pathways in a more gen¬ 
eral open quantum framework. Often one may want to 
select a particular pathway without also addressing un¬ 
wanted transitions. When the spin transitions are close 
or there is an overlap between the desired excitation fre¬ 
quencies and other transitions in the system, sophisti¬ 
cated control techniques are required. 

In this work we applied optimal control theory (OCT) 
to find control sequences that select one DNP pathway 
to produce the largest possible nuclear polarization en¬ 
hancement. OCT is a valuable tool in magnetic reso¬ 


nance and for control of quantum informatiorP^. Opti¬ 
mal control pulses have been used previously to demon¬ 
strate universal control in electron-nuclear systems with 
microwave irradiation onlj®l. 

We use DNP to demonstrate the improved control 
available in open quantum systems through optimal con¬ 
trol methods. For unitary systems, there is always a 
perfect solution with optimal control, even in systems 
with a distribution of Hamiltonians such as powder sam¬ 
ples or samples in the presence of large field inhomo¬ 
geneities. Finding optimal solutions becomes more com¬ 
plex when decoherent processes are involved. There 
has been theoretical work applying OCT to open quan¬ 
tum systemPilti^, including optimization of DNP pro¬ 
cesses with simultaneous RF and microwave controP^HSl^ 
Others have also studied optimal control in the con¬ 
text of quantum measurementfi^, quantum systems cou¬ 
pled tonon-Markovian environmental, and minimizing 
nois^miSl Nonunitary c ontro l of DNP has also been con¬ 
sidered in quantum dotP^l^. We present here an exper¬ 
imentally realized a scheme for microwave only optimal 
control of DNP in an electron-nuclear system with relax¬ 
ation mechanisms included. 

Our system consists of one electron hyperfine coupled 
to one hydrogen nucleus, forming a four level system (Fig. 
[^. There are two available DNP pathways, depending on 
which transitions are excited. Applying microwave irra¬ 
diation on resonance with the zero quantum transition in¬ 
duces an electron-nuclear mutual spin flip-flop, which di¬ 
rectly transfers polarization from the electron spin to the 
nuclear spirP. This is known as the solid effect (SE). Al¬ 
ternatively, the Overhauser effect (OE) indirectly causes 
this flip-flop through a cross relaxation process that oc¬ 
curs when both electron resonances are saturated (see 
Fig. 0 for diagrams of both mechanismsThe zero 
quantum transition is ordinarily forbidden, and cross re¬ 
laxation {Tx in Eig. is not allowed. Both the solid ef¬ 
fect and Overhauser effect can occur, however, when the 
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tonian, which contains a hyperfine coupling with both 
isotropic (^) and anisotropic (B) terms. In the high 
field limit the electron Zeeman term will be much larger 
than the other terms, we take the secular approximation 
with respect to the electron spin quantum number and 
the drift Hamiltonian is 


(a) 


a|H>+P|H> (b) a|H>+PIH> 


Hdrift — <-^sSz + W/Jjz + ASzIz + BSzIx (1) 


FIG. 1. The energy level diagram on the left illustrates the 
Overhauser effect: the two electron resonances, r'i and v^, 
are saturated, and polarization is transferred through the zero 
quantum cross relaxation, T^. The diagram on the right rep¬ 
resents the solid effect: the zero quantum transition, is 
directly irradiated. 

nuclear Zeeman states are mixed due to an anisotropic 
hyperfine coupling. 

Under continuous wave (CW) irradiation all transi¬ 
tions are excited to varying degrees, driving both pro¬ 
cesses. The OE leads to a positive nuclear polarization, 
and the SE to negative nuclear polarization. Conse¬ 
quently if both the OE and SE are driven, the net polar¬ 
ization enhancement is less than the polarization achiev¬ 
able by either process acting on the system alone. It is 
therefore advantageous to find a control sequence that se¬ 
lects one DNP process and suppresses the other. A pulse 
sequence capable of driving exclusively one of these two 
opposing methods of DNP must account for both deco- 
herent processes and unitary evolution. 

The theoretical nuclear polarization enhancement is 
limited by the ratio 7e/7n ■ For hydrogen, the max¬ 

imum achievable enhancement is 660 . In any physical 
system the enhancement will be less than the theoretical 
limit due to a factors such as imperfect saturation, leak¬ 
age effects from relaxation, and the competition between 
various DNP processes. 


II. QUANTUM DESCRIPTION OF DNP 


where S and / are the spin operators for the electron 
and the nucleus respectively, ujs/uji the electron/nucleus 
Larmor frequency. For efficient polarization transfer, it 
is necessary to use a system in which the anisotropic hy¬ 
perfine is on the same order as the nuclear Zeeman term, 
so that the eigenstates are mixtures of nuclear spin up 
and spin down and the zero quantum transition is weakly 
allowecP. 

The control Hamiltonian for the OE is produced by 
microwave irradiation on resonance with the electron spin 
flip transitions. In the rotating frame of the electron this 
is: 

Hoe = ^dSx O 1 (2) 

where ujd is the Rabi frequency of the drive field. Simi¬ 
larly for the SE, the microwave control on resonance with 
the zero quantum transition yields the control Hamilto¬ 
nian: 

HsE = UJd (Sxlx + Syly) ( 3 ) 

The remaining operators needed to mathematically de¬ 
scribe DNP are the non-unitary operators that drive the 
decoherent processes. The SE requires only Tf relax¬ 
ation, while both Tf and Tx are necessary for the OE. 
All the unitary and non-unitary operators for the solid 
and Overhauser effects are listed in Table HI 

To illustrate that these operators are sufficient, we con¬ 
sider the following discrete maps which increase the nu¬ 
clear polarization. Here we assume that Ti of the nuclear 
spin is much longer than Tf and can be neglected. 

For the Overhauser effect: 


Historically DNP has been described using a semi- 
classical model with systems of rate equations. We have 
the tools, however, to analyze the polarization transfer 
as a fully quantum process with system-environment in¬ 
teractions. We can write both the solid effect and Over¬ 
hauser effect as quantum maps, and we do so here using a 
Kraus operator representation. The minimum set of op¬ 
erators necessary to determine these maps are the drift 
Hamiltonian, Hdrift, the control Hamiltonian, Hcontrol, 
a set of Kraus operators for the T[ relaxation of the elec¬ 
tron, and Kraus operators for Tx, the zero quantum cross 
relaxation between the electron and the nucleus. 

The drift Hamiltonian is the same for both DNP mech¬ 
anisms. Hdrift is the sum of the Zeeman interactions for 
the electron and the nucleus and the hyperfine coupling 
between the two spins. We assign a particular Hamil- 


Aoe = [Uoe] ^ [Tx] ^ [Tf] (4) 

1 - Uoe saturates the electron spins, removing any 
electron magnetization. 

2 . Tx; drives population to the nuclear spin |4,) state. 

3 . Tf removes entropy from the system. 

And the Solid effect: 

Ase = [Use] -t [T^] (5) 

1 - Use transfers polarization from the electron 
through the zero quantum transition. 

2 . Tf removes entropy from the system and resets the 
electron spin. 
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The other set of Kraus operators that is required to 
describe DNP is that which implements the electron- 
nuclear cross relaxation, (see Fig. [^. This is es¬ 
sentially a Ti process that acts on the zero quantum sub¬ 
space, a depolarizing channel that leads to transitions 
II j5) |t a-*-). As in the electron Ti case the system 
should return to thermal equilibrium. The following op¬ 
erators satisfy these requirements: 


TABLE I. Maps for both mechanisms in rotating frame of the 
electron. 


In a physical system decoherence is continuous, and 
we cannot sequentially apply the unitary processes then 
decoherent processes compose these discrete maps. In¬ 
stead we can move to a continuous time map and find 
the generators for DNP. We choose the Kraus operator 
formalism to describe the decoherent processes and de¬ 
termine the continuous map. 

In the Kraus operator form, a set of operators {Mk} 
act on the state p yielding the map, 

( 6 ) 

k 

This map is completely positive and trace-preserving if 
the Kraus operators satisfy the condition, 

J2mIm,=1 (7) 

k 

We start with the Kraus operators that describe Ti re¬ 
laxation. During the Ti process, the electron spin flips 
while the nuclear spin remains in the same state (indi¬ 
cated by Tf Fig. [^. A set of Kraus operators driving 
this type of relaxation must act as the operator 0 1/ 
in the energy eigenbasis. We also require the Kraus op¬ 
erators for a Ti process to return the system (or relevant 
subsystem) to thermal equilibrium. The following set of 
Kraus operators describes Ti relaxation on the electron 
spin: 
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( 9 ) 


with Pzq = exp{—h(ujs — uji)/ 2kBT) and Czq = 
exp(—t/T^g), where Tzq is the zero quantum cross re¬ 
laxation. These Kraus operators satisfy the condition 

We now have the Kraus map for the total evolution 


P = 








( 10 ) 


which we rewrite in terms of a single set of Kraus oper¬ 
ators, {Mk}, 


p' = Y,MkpMl ( 11 ) 

k 


.4, = Vp(” '^)®I 

A 3 = \/l -p ^ ^ 

( 8 ) 

where the parameter e = exp(—t/Tf) sets the rate of the 
process and p = exp{—hu!s/2kBT) determines the hnal 
polarizatioipS. 


The next step is to find the reduced map that acts on 
the nuclear spin subsystem and to write it in the Kraus 
form. We start by transforming the Kraus map to a 
supermatrix representation that acts on the vectorized 
density matrix, p, 

P = S3 (12) 

In a general sense, we can consider the electron to be the 
environment of the nucleus. We evaluate the reduced 
map on the nuclear spin by tracing over this “environ¬ 
ment” . This is accomplished by applying S to a known 
state of the environment and taking the partial trace, 

Sn = Trsi^PE ® Pn) 


(13) 
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The Kraus operators acting on the reduced space can 
be found from the Choi matrix, which is given by a 
reshuffling of the reduced supermatrix, Sj^. The Choi 
matrix is defined as 

Ac = (I O A) ^ ® Eij = E,j (g) A(Eij) (14) 

ij ij 


where Eij = \i) {j\ and the vectors {|*), |j)} form a basis 
of the Hilbert space. The procedure for extracting Kraus 
operators from the Choi matrix follows. We rewrite Ac 
in terms of its normalized eigenvectors, \ak), 


Hamiltonian parameters and Ti as follows, 

hujs 


r = 7 i = exp 


2kBT 


(17) 


a = 


(1 — e */^i) 


A+ = 


, (T4-3y2)(t-4Ti) 

P± - , -^ + ^(t ) 

8V3 ± 2y2ri 

—STi + t (2 i + i4Ti (7y_ + 7 +) T 2 '\/ 27 i) 

" ^ 2(t-4ri) 

+ 0 (^ 2 ) 


with 


Ac = ^ \ak) {ak\ (15) 

k 


The columns of Ac form a set of columnized Kraus op¬ 
erators {Ak} satisfying the condition of Eq. 0 with 
d? operators in the set, where d is the dimension of the 
Hilbert space. 

Using this procedure we find the following Kraus op¬ 
erators for the nuclear spin map. The final Kraus oper¬ 
ators for both DNP processes are similar in form to the 
Ti Kraus operators of Eq. Q and do in fact describe po¬ 
larizing channels on the nuclear spin. For the solid effect 
effect the four Kraus operators found have the structure: 


' Vo 0 y 


(16) 


V± 


+ 4B'^ ± AAuJi + uj 


(18) 


/3± and A± are approximated to first order in time, t. 
This approximation is valid under the assumption that 
the time step is much shorter than the relaxation time, 
Ti. 

Similarly for the Overhauser Effect, we find the re¬ 
duced Kraus operators. 




,/o v/ri(i-r.)^ 

Vo 0 J 
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0 0\ 
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AV 0\ 

0 0 


(19) 


with the parameters, 


ri = 1-271 (l-e‘/^0 

h{ujs - UJi) 

r.=7. = exp 

' ^ g-qi/Ti + l/T,,) 


( 20 ) 


where the parameters T, a, /3±, and A± depend on the 


and 


P± = 


'-t T 4T, ± ^4 X (-4tT, + 8 ry + + 2T^ir,_ + ry+)2 - 27,(1 - 7 ,))) 

4T, 


2 i 


A'± = 


i (±t(l - 27 ,) + v'4 X (-4tT, + 8 ry + t2(i + 2T^{rj_ + r,+ )^ - 27 , + 27 ) 

— t{2i + 4:Tx{ri- + rj^)) 


0{t^) 

+ o{e) 


with 7 i and ry as defined in Eqs. 0 and (|18[) respec- tively. Again, for simplicity, we have taken and A(j_ 
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FIG. 3. Timing diagram of the experiment: the n/2 satura¬ 
tion is applied for a variable length of time followed by the 
NMR readout pulse. 


FIG. 2. The four possible transitions in the two spin system. 
1 and 2 are electron only transitions. 3 is the zero quantum 
transition, and 4 the double quantum. 


to first order in t with the assumption that t is shorter 
than Ti and T^. If we make the same approximation for 
Fi and F^:, then the cross relaxation term dominates and 
the Kraus operators M[ and take the same form as 
those for the solid effect {Mi and M 2 in Eq. but 

with the direction of polarization reversed: 


M[ 

M2 —>■ 



0 

00) 

0 0 \ 

Vm 0 ) 


( 21 ) 


These Kraus operators satisfy the requirement of Eq. Q 
for a CPTP map. {Mk} and {M'f.} resemble the Kraus 
operators for a return to thermal equilibrium, each with 
a new equilibrium defined by 71 and ■jx- These reduced 
maps on the nuclear spin drive the spin to a hyperpo- 
larized state, and with opposite sign for the solid effect 
versus the Overhauser effect. 


III. OPTIMAL CONTROL THEORY 

The generators given above assume perfect control. In 
reality, we cannot produce the perfect unitaries Use and 
UoE- Any microwave field that we apply to the system 
will excite all four transitions to some degree, resulting 
in non-ideal unitaries that will drive the Overhauser ef¬ 
fect and solid effect simultaneously and reduce the final 
nuclear polarization. In order to increase selectivity, we 
incorporated optimal control theory (OCT). 

OCT allows the design of pulse sequences that per¬ 
form any desired unitary on the four level system. There 
are many unitary operators that result in a hyperpolar- 
ized nuclear state. Exciting the transitions that drive the 
Overhauser effect, for example, requires a unitary that 
acts only on the subspace of the electron spin, whereas 
the solid effect is driven by unitary operators acting on 
the zero quantum subspace. Eor the experimental im¬ 
plementation we choose a single operator to optimize to¬ 


wards; in principle we can optimize to any operator on 
the electron subspace that leads to a polarization trans¬ 
fer. 

There are four possible transitions in the hyperfine cou¬ 
pled electron-nuclear system, as shown in Eig. De¬ 
pending on which transitions are driven during the DNP 
experiment, the effective unitary will lead to a combina¬ 
tion of the OE and SE. As will be explained in the next 
section, the pulse sequence used to drive DNP in this 
work is a saturation train (Fig. [^, where each pulse will 
either be a hard pulse or a composite OCT pulse. The 
angle of rotation that a single pulse performs on each 
transition in Fig. determines the degree of saturation 
of that transition. 

By varying the rotation angle of each transition we can 
explore the space of DNP and find points in this space de¬ 
termined that yield a nuclear polarization enhancement. 
Fig. [^presents the DNP enhancement of the nuclear spin 
signal for varying saturation levels of the four transitions. 
There are four angles 9i corresponding to the transitions 
of Fig. These plots illustrate several representative 
slices of the fnll space of DNP. The x-axes of (a)-(c) in¬ 
dicate the rotation angle of the pulses in Fig. [^on tran¬ 
sitions 1 and 2 (i.e. the transitions that involve flipping 
only the electron), with no excitation of the zero quan¬ 
tum and double quantum transitions. These points in 
the space indicate unitaries which drive the Overhauser 
effect only and lead to positive nuclear polarization en¬ 
hancement. When the cross transitions (3 and 4 in Fig. 

are also excited (increasing saturation corresponds to 
larger rotations on the y-axes), the microwave control 
also induces the solid effect and drives the nuclear spin 
towards a negative polarization. In regions where sin¬ 
gle quantum transitions and cross transitions are both 
excited, the net nuclear polarization is less than the the¬ 
oretical maximum of 660. 

Experimentally it is difficult to drive polarization en¬ 
hancement through the solid effect. The zero quantum 
and double quantum transitions are only weakly allowed, 
and therefore require large microwave powers to saturate. 
In this work we choose instead to drive the Overhauser 
effect by optimizing tt/2 rotations on the 1 and 2 tran¬ 
sitions with no excitation on the cross transitions, 3 and 
4. This case is indicated by arrows on figures (a) and (c) 
of Fig. ^ 
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FIG. 4. Color map indicates the DNP enhancement, with the maximum of 660 in red and minimum of -660 in blue. The x 
and y axes indicate the degree of rotation for a repeated train of pulses (as in Fig. i) on the transition given by the axes 
labels. The angles 6i correspond to the transitions in Fig. The following cases are depicted: (a) Increasing saturation on 
transitions 1 and 3 while applying 7r/2 pulses to 2; (b) Increasing the 2 and 3, with 7r/2 pulses on the other electron resonance, 
1; (c) increasing 2 and 4 with full 7r/2 saturation on 1; (d) increasing 3 and 4 with both electron resonances saturated; and 
(e) increasing the electron resonance simultaneously and increasing 4 with 7r/2 pulses on 3. The ideal case for the Overhauser 
effect is occurs at the points indicated by arrows in (a) and (c). 


IV. PULSE FINDING 


We accomplish the saturation of the electron reso¬ 
nances for the OE through a 7r/2 microwave pulse train 
as shown in Figure The goal of the pulse finder is to 
produce composite pulses that perform the 7r/2 rotation 
on both electron resonances without any rotation on the 


zero quantum (or double quantum) subspace. The de¬ 
sired unitary for Overhauser DNP is 

U = exp(—f7r/2S'a; ® 1) (22) 

Our pulse optimization method uses the Nelder-Mead 
simplex algorithm to numerically search over the control 
parameter spac^^. The pulses are restricted to on/off 
modulation; the total pulse consists of several time steps 
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during which the microwave control is either on or off 
and the lengths of each step are left as parameters. 

To optimize the pulses, we first find the map that the 
OCT pulse performs when relaxation is present in the 
system and compare this to the map produced by the 
desired 7r/2 rotation. We have a set of operators {Mk} 
such as in Eq. Q for each time step. The successive 
application of these Kraus operators for each segment of 
the pulse gives the map for the composite pulse: 

P'=E< ■■■ 

ni L 

In order to speed up the computation time, we use the 
procedure from Section |III| to find the smallest set of 
Kraus operators from the Choi matrix. 

The pulse finder’s fit function measures the overlap of 
the superoperator generated by the OCT pulse wit h the 
desired unitary, represented by the gate fidelitjEHH^ 

F = ^ Tr{U^Mk)^/2‘^ (24) 

k 

where U is the unitary we want to perform and {Mk\ is 
the non-unique set of Kraus operators that describe the 
pulse superoperatoi!^. 

We found two sets of OCT pulses: one including the 
electron Ti and the cross relaxation processes (Tj,), and 
one considering only unitary evolution. The pulse found 
with optimal control but without including relaxation 
consists of two pulses (p„) with a delay (r„) between 
them: 

(n) - (pi) - (t 2) - (P2) - (rs) (25) 

This composite pulse, which was found neglecting relax¬ 
ation, has a fidelity with the desired unitary of 0.56. 

The OCT pulse found by including relaxation pro¬ 
cesses is three square pulses with delays as follows: 

(n) - (Pi) - (t2) - (P2) - (ts) - (P4) - ( 7 - 4 ) (26) 

The computed fidelity for this pulse is 0.74. It appears 
that the additional square pulses in the OCT sequences 
correct errors that hard pulses do not. These fidelities 
may appear low, but it should be noted that they cor¬ 
respond to a single application of the composite pulse. 
We need only to saturate a transition in order to trans¬ 
fer polarization. Therefore the angle of the pulse does 
not matter as it is repeated many times in the saturation 
train, and we can still achieve large polarization enhance¬ 
ments. Additionally, these fidelities are calculated for the 
unitary that acts on the full electron-nuclear space. For 
DNP we are not concerned with the full map that acts on 
the electron-nuclear system, but we are interested only in 
the action of the control pulses on the nuclear subspace. 
Over the length of the saturation train the reduced map 
on the nuclear spin after tracing out the electron is essen¬ 
tially a polarizing channel for all sets of the pulses found. 



OE 

SE 

Open OCT 

Closed OCT 

Hard Pulses 

II 

0.25 

0.25 

0.25 

0.25 

0.25 

IX 

0 

0 

8.13e-8 

-9.25e-8 

5.46e-8 

lY 

0 

0 

-1.43e-7 

-1.76e-7 

5.53e-ll 

IZ 

1.88e-4 

-1.85e-4 

1.79e-4 

1.24e-4 

5.90e-5 

XI 

0 

0 

-1.69e-5 

-3.98e-5 

3.88e-5 

XX 

0 

0 

-6.29e-12 

1.38e-ll 

8.80e-12 

XY 

0 

-1.28e-6 

9.93e-12 

2.94e-ll 

1.36e-14 

XZ 

0 

0 

-1.25e-8 

-2.06e-8 

9.61e-9 

YI 

5.13e-6 

0 

2.64e-6 

2.68e-6 

-5.70e-6 

YX 

0 

1.28e-6 

5.Ole-13 

-8.52e-13 

-1.07e-12 

YY 

0 

0 

-1.02e-12 

-2.04e-12 

-3.43e-14 

YZ 

3.85e-9 

0 

2.07e-9 

1.44e-9 

-1.27e-9 

ZI 

-6.52e-6 

-1.88e-4 

-1.35e-5 

-6.66e-5 

-1.30e-4 

ZX 

0 

0 

-4.29e-12 

2.51e-ll 

-2.94e-ll 

ZY 

0 

0 

9.Ole-12 

4.92e-ll 

6.99e-15 

ZZ 

1.20e-5 

1.40e-7 

-9.92e-9 

-3.42e-8 

-3.16e-8 


TABLE II. Final states given by the superoperators for the 
Overhauser effect and solid effect compared to the states 
driven by the superoperators for the three sets of pulses found. 

This map takes the identity state to a polarized state as 
follows, 

K{t) = t+pZ (27) 

with the value of p determining the final nuclear polariza¬ 
tion for each type of pulse. The OCT pulses found includ¬ 
ing the full open system dynamics act as the strongest 
polarizing channel and have a fidelity of 0.95 with the re¬ 
duced map produced by the ideal Overhauser effect. In 
contrast the fidelity of the reduced system dynamics with 
the ideal case is 0.66 for the closed system OCT pulses 
and 0.31 for hard pulses. 

To analyze the map on the full space for each of the 
pulses, we rewrite it in the supermatrix formalism. Fig. 

gives the superoperator matrices, S, in the Pauli ba¬ 
sis for the solid effect, Overhauser effect, both sets of 
OCT pulses, and hard pulses. The eigenvector of S with 
eigenvector 1 is equal to the final state into which the 
superoperator drives the system. These final states for 
the Overhauser effect and solid effect superoperators are 
given in Table [TT| The IZ and ZI terms are proportional 
to the final nuclear and electron polarizations, respec¬ 
tively. 

The final state coefficient indicate that the OF drives 
the system to a positive nuclear polarization while depo¬ 
larizing the electron spin. In contrast, the SE results in a 
negatively polarized nuclear state but does not perform 
a strong depolarizing channel on the electron spin. Ad¬ 
ditionally the SE creates coherences in the zero quantum 
terms, XY and YX, while the OF only leads to single 
quantum coherences in the final density matrix. 

From the superoperator matrices in Fig. [^we can see 
qualitatively that the superoperator for OCT pulses has 
the same structure as the superoperator for the Over¬ 
hauser effect with perfect 7r/2 rotations on the electron 
only transitions. The closed system OCT pulses and hard 


(23) 














pulses are also similar in structure to the OE map, but 
have increasing contributions from terms not present in 
the OE map. 


V. SIMULATIONS 

Figure shows the results of simulating the OCT and 
non-OCT pulses. The simulations show that both sets 
of OCT pulses (solid black for open system OCT and 
dot-dashed red for closed system OCT) give a greater 
enhancement than hard 7r/2 pulses (dashed blue). The 
open system OCT pulses yield a nuclear polarization of 
654, almost the full enhancement. The closed system 
OCT pulses produce a nuclear polarization 491, while 
the enhancement due to hard pulses is 140. 

We simulated the performance of the open system 
OCT pulse with variations in both Rabi frequency and 
the anisotropic hyperfine coupling and compared it to 
hard pulses with the same variations (see Fig. [^. The 
OCT pulse was optimized for a Rabi frequency of 8 MHz, 
the frequency experimentally achievable at a reasonable 
power. The simulations show that the OCT pulses are 
quite robust to changes in Rabi frequency, while hard 
pulses reach the maximum polarization enhancement at 
a Rabi frequency equal to anisotropic hyperfine coupling 
(14 MHz in this simulation). The enhancement falls of 
relatively quickly for hard pulses as compared to the OCT 
pulses. 

Fig. [3b) presents the results of varying the anisotropic 
hyperfine frequency, B (as in the Hamiltonian of Eq. 

The OCT pulses give the highest enhancements when 
the hyperfine coupling is near maximal, where there is 
the largest amount of nuclear state mixing. The hard 
pulses actually improve with lower anisotropic hyperfine 
values. This result is somewhat nonphysical, however, 
as the cross relaxation rate is the same for all values of 
B. In reality the cross relaxation is weaker when B is 
smaller. In these simulations where is constant and 
B decreases, the hard pulses improve because the solid 
effect becomes less likely as the zero quantum transitions 
is more strongly forbidden, while the Overhauser effect 
is still allowed through the constant cross relaxation. 

One possible leakage pathway in the physical system is 
the double quantum transition. While relaxation through 
flip-flip interactions will be a weaker process than the flip- 
flop that drives DNP, both pathways are present. We 
have simulated DNP in a system with a double quantum 
relaxation rate of 0.5 times the zero quantum rate. The 
results are shown in Fig[^ When double quantum relax¬ 
ation is included, the DNP enhancements are lowered for 
all pulses: 218 for open OCT, 152 for close OCT, and 38 
for hard pulses. 

It should be noted that there is a significant differ¬ 
ence between the simulated system and the physical sys¬ 
tem. In the simulations we are looking only at the two 
spin system consisting of one electron and one proton 
with the couplings found in malonic acid. We find the 


DNP enhancement of the coupled proton. In the experi¬ 
ment, however, we observe the bulk proton signal, as the 
number of protons coupled directly to an electron is too 
small to detect directly. The actual experimental nuclear 
polarization depends on the couplings between the bulk 
protons and the hyperfine coupled protons, and the build 
up times are limited by spin diffusiorP^. 

VI. EXPERIMENT 

The system used in this work is irradiated malonic 
acid, an organic radical with one electron hyperfine cou¬ 
pled to one proton. Malon ic aci d is a well characterized 
sample widely used in ESRP^Iin! The malonic acid sam¬ 
ple used was a single crystal irradiated for 5 hours with 
8 keV X-rays and subsequently annealed at 45° C for 15 
hours. 

The Hamiltonian for this two spin system is the drift 
Hamiltonian of Eq. Q. In an external field of 3406 G, 
the electron Larmor fr equen cy is 9.59 GHz, the nuclear 
resonance is 14.57 MH #° l ^d The terms A and B depend 
on the orientation of the crystal in the external magnetic 
field. In the orientation that maximizes the mixing of 
the nuclear states, the isotropic (A) and anisotropic {B) 
parts of the hyperfine coupling are -42.7 MHz and 14.7 
MHz respectively. We determined this orientation in a 
GW-ESR spectrometer. 

We built a double frequency probehead with an ESR 
resonator surrounded by a split coil for NMR detection. 
The schematic for the probe is shown in Figj^ The ESR 
resonator is based on the standard loop-gap resonator de¬ 
sign, but is easier to construct than typical bridged loop- 
gap resonators commonly used in DNP experimentP^Ml 

The resonator is soldered onto the outer conductor of 
a microwave coax at one end, while the other end is open 
(see top view of Figj^. A tuning screw on the opposite 
side of the shield allows tuning and matching of the ESR 
resonance over a 1 GHz range around 10 GHz. 

We performed three sets of experiments, all using the 
basic pulse sequence given in Figure with hard pulses, 
closed system OCT, and open system OCT pulses as 
the 7r/2 saturation pulse. To measure the thermal sig¬ 
nal without DNP, we used an Ernst angle detection with 
400,000 scans. Due to the large dipolar coupling between 
protons in malonic acid, T 2 of the proton signal was less 
than the deadtime of the probe. The NMR signal was 
detected with a magic echo sequence to overcome this 
issue. The proton linewidth was 40 kHz. 

VII. RESULTS 

Figure [T0| shows the build up curves for the three sets 
of pulses. The open system OCT pulses produce the 
highest nuclear polarization, 11.1 times the enhancement 
produced by hard pulses, and 2.2 times that produced by 
the OCT pulses found without relaxation (closed system 
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FIG. 5. Supermatrices for OE (a) and SE (b) as well as for the three sets of pulses: open OCT (c), closed OCT (d), and hard 
pulses (e). Dark Red indicates elements with magnitude 1, and at the other side of the spectrum dark blue elements are equal 
to 0. The operator describing the open OCT pulses has the best overlap with that for the OE. The two other pulses have large 
contributions from components not present in (a), indicating that they are less effective at isolating the OE 
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FIG. 6. Comparison of OCT pulses to hard pulses. Black 
solid line is saturation by open system OCT pulses, dashed 
red is closed system OCT pulses, and short dashed blue is 
saturation with hard pulses. 


pulses). The closed system pulses yield a polarization 
5.1 times that of hard pulses. These ratios are consis¬ 
tent with those found in the simulations. The maximum 
polarization enhancement achieved with the open OCT 
pulses was a factor of 180 ± 36. 

From fits of the data in Figure 10 we find the build 
up times for each curve. For all three curves the rate 
determining process is spin diffusion, as we detect the 
bulk proton signal and do not measure the polarization 
of the spins directly coupled to the electron defects. The 
hard 7r/2 pulses produce a build up time of 12.7 ± 1.8s , 
the closed OCT build up is 14.9 ± 1.7s, and 13.2 ± 1.5s 
for the open OCT pulses. 

Figure El shows the measured effective T) of the bulk 
nuclear spins after DNP. A fit of this curve to a single 
exponential decay gives Ti = 10.5 ± 1.6 s. Like the DNP 
process, the return to thermal equilibrium depends on 
spin diffusion. The fast electron Ti process will depolar¬ 
ize nuclear spins in its vicinity, which will then diffuse 
through the bulk. Thus we expect the effective Ti of the 
bulk nuclear spins to be on the same order as the polar¬ 
ization build up times. We attribute the slight difference 
between the polarization and depolarization times to the 
fact that there are more pathways for the bulk spins to 
relax through than there are pathways for polarizing. 

These findings provide an example of how optimal con¬ 
trol theory can improve polarization transfer in DNP ex¬ 
periments, and, more generally, provide superior control 
in open quantum systems. Even with a relatively sim¬ 
ple pulse sequence with few time steps we can achieve an 
order of magnitude increase in polarization using OCT 
pulses versus hard pulses. 

While the OCT pulses do produce a nuclear polariza¬ 
tion enhancement of 180, this increase is less than the 
theoretical maximum enhancement of 660 and less than 
the simulated value. The inability to reach this level of 
enhancement could be due to several factors already dis- 


(a) 




Anisotropic Hyperfine (MHz) 

FIG. 7. Simulations of pulses varying the Rabi frequency (a) 
and the anisotropic hyperfine coupling (b). The solid black 
line indicates simulations of the OCT pulses, red dashed is 
hard pulses. 


cussed in Section |Vj such as leakage pathways or asym¬ 
metric saturation of the electron resonances. The simula¬ 
tions presented earlier have shown that leakage through 
double quantum relaxation (in Fig. |^, or asymmetric 
saturation (|^ )will also reduce the final nuclear polariza¬ 
tion. We have also seen that the hard pulses are more 
susceptible to such deviations from the idealized case, 
and it is not surprising that the hard pulses are worse 
than the simulations by a factor of 10, while OCT pulses 
are only a factor of 5 different than their simulated value. 


VIII. DISCUSSION 

These results could be extended to find control se¬ 
quences for larger spin systems. In particular we have 
neglected any nuclei further from the electron defect than 
the nearest proton. These nuclear spins may still have 
small hyperfine interactions with the electron spin, and 
we have not allowed any transfer pathways in the pulse 
optimization which might polarize these spins directly. 
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FIG. 8. DNP with double quantum relaxation. The black 
solid line is open OCT pulses, red dashed closed OCT, and 
blue dotted is hard pulses. Here Tdq = 2Tzq. 


FIG. 11. Effective Ti measurement after polarization for 10s. 
Fit to a single exponential given by black line. Ti measured 
to be 10.5 ± 1.6s 


ESR resonator 




FIG. 9. Probe Diagram: the figure on the left shows the side 
view of the probe, and on the right is the top view. The ESR 
resonator in the shape of a loop-gap is concentric to the NMR 
coil. 



Time (s) 


It may also be useful to design control pulses for sam¬ 
ples other than single crystals. In the case presented here 
the two spin system was assumed to be in the same ori¬ 
entation throughout the sample. It may be particularly 
useful to consider DNP in powder samples. This would 
require optimizing pulses for a set of internal Hamiltoni¬ 
ans. 

This method of pulse finding is useful for quantum con¬ 
trol in systems where relaxation plays a strong role in the 
dynamics. The Kraus operator description is particu¬ 
larly convenient as one can write a modular pulse finding 
code in which decoherent processes can be added to the 
superoperator successively. In principle another super¬ 
operator formalism, such as the Lindblad representation 
under the condition that the system is Markovian, should 
give the same results. Including relaxation processes in 
the pulse finder produces pulses that select the desired 
transitions more accurately than using a superoperator 
that accounts only for unitary processes. This experi¬ 
ment demonstrates that we can design control sequences 
for open quantum systems that account for interactions 
with the environment. 

Not only is it possible to create control operations in 
dissipative systems, but we have shown that an open 
quantum system approach can provide significant im¬ 
provements over unitary control. We anticipate that su¬ 
peroperator optimized pulse design will lead to superior 
control in spin based systems and other quantum devices. 


FIG. 10. Results of the experiment with fits of the buildup 
curves. Black triangles are the data from saturation with open 
system OGT pulses, red circles are for closed OGT pulses, and 
blue squares are for hard pulses. 

We could also include the dipolar couplings between nu¬ 
clear spins to simulate spin diffusion, the rate limiting 
step in this experiment. 
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